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Properties of some nonlinear Schrodinger equations 
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Abstract. We update our understanding of nonlinear Schrodinger equations motivated 
through information theory. In particular we show that a q— deformation of the basic nonlinear 
equation leads to a perturbative increase in the energy of a system, thus favouring the simplest 
q = 1 case. Furthermore the energy minimisation criterion is shown to be equivalent, at 
leading order, to an uncertainty maximisation argument. The special value rj = 1/4 for the 
interpolation parameter, where leading order energy shifts vanish, implies the preservation of 
existing supersymmetry in nonlinearised supersymmetric quantum mechanics. Physically, i] 
might be encoding relativistic effects. 
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1. Introduction 

We have previously argued, see for example [T], that the maximum uncertainty (entropy) method 
[2], used widely in several disciplines is also a plausible framework to investigate potential 
nonlinear generalisations of Schrodinger's equation [4j. 

We consider here for simplicity a single particle in one space dimension. The basic nonlinear 
equation of Ref . [1] is 
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p{x) = 'ilj*{x)'tp{x) and p±{x) = p{x it rjL). The nonlinearity length scale L and the energy 
parameter £ are constrained through the relation £L^ = fi^/Am. The dimensionless parameter 
< rj < 1 plays the role of an interpolating parameter and also regularises potential singularities. 

The nonlinear equation preserves a number of important properties of the linear Schrodinger 
equation, such as the conservation of probability, existence of the usual plane wave solutions 
and invariance under a scaling of the wavefunction, ^p — > Xip. 



2. Physical meaning of the r/ parameter 

In Ref. [5] nonlinear Dirac equations were constructed using an axiomatic approach [Gj . In taking 
the non-relativistic limit of those equations [7j it was found that the potential singularities in 
the naive limit were regularised if one kept sub-leading relativistic corrections coming from the 
lower components of the Dirac spinor. This suggests an interpretation in which the rj parameter 
in the nonlinear Schrodinger equation ([T]) is encoding relativistic effects present already in the 
linear quantum theory, while the nonlinearity is parametrised by the scale L. 

Alternatively, one may proceed as follows. In Ref.[8j we computed the perturbative effect of 
the nonlinearity ([T]) on the energy spectrum of the linear theory with smooth external potentials 
V{x). For an unperturbed state with quantum number n, with nodes in the wavefunction, 
the result is 



where a is the characteristic length scale of the linear theory. The coefhcients Cnp depend on the 
slope of the unperturbed wavefunctions near the nodes, and hence on the form of the potential. 
Notice however that the dependence of the energy shift on rj is universal, that is, independent 
of V{x). (For states without nodes the energy shift is of lower order, 0{L/a)^). 

It was suggested in |T| that rj be fixed by minimising the energy of the system. The result ([3]) 
indicates that this happens when -q ~ 0.80. Now, to leading order, the energy shift ^ is given 
by < F{p) >, evaluated using the unperturbed states. Assuming integration over an infinite 
domain, this reduces after a simple shift of variables tij to Iatl — /p, where Ip is the Fisher 
information measure that leads to the linear Schrodinger equation Oin] and Inl = —Ikl is the 
information measure that leads to the nonlinear generalisation ([T]); Ikl is the Kullback-Liebler 
measure. 

Thus minimising the energy shift is equivalent to minimising the information content of the 
nonlinear theory relative to that of the linear theory, and this is achieved for excited states at 
the universal value rj ~ 0.8. 

An interesting aspect of the result ([3]) is that it vanishes for values r/ = 0, 1/4, 1. The r/ = 
case is obvious becuse then the nonlinearity is switched off. The ry = 1 case is amusing because 
there the nonlinearity exists but the intermediate perturbative calculations are ill-defined if one 
starts with r/ = 1; thus r/ ^ 1 may only to taken after the calculations. We will see in the next 
section that the vanishing of the energy shift at ?7 = 1 is likely to be an accident, valid only at 
leading order. 

So what about the ry = 1/4 case? This is quite puzzling and we still do not understand if 
it indicates some hidden symmetry of the theory at that value. However we see that at that 
special value the nonlinearity can help preserve, to leading order, an existing supersymmetry 
of the linear theory: Consider one dimensional super symmetric quantum mechanics [TT] with 
two partner Hamiltonians Hi , H2 ■ Except for the ground state of Hi the two Hamiltonians are 
isospectral but with differing wavefunctions. For example, the first excited state of Hi will have 
the same energy as the ground state of H2; the corresponding wavefunction of Hi will have one 
node while that for H2 will be nodeless. Thus if we nonlinearise both Hamiltonians as in ([T]), 
we see that at rj = 1/4, all energy shifts, whether for excited or ground states will be of order 
{L/a)^ and so to leading order the pair of Hamiltonians will have their original supersymmetry 
preserved. 

3. Robustness of results 

The nonlinearity in Ref. [4J was obtained using, arguably, the simplest generalisation of the Fisher 
information measure which preserved various desirable properties of the former. How sensitive 
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are the above results to the particular choice of /xl? We can study this issue by considering 
deformations of the Kullback-Liebler measure. 

Consider therefore the g— deformed measure {q > 0), 



M(ri, q) = — J / dx dt pix) 
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with 
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This measure reduces, when g = 1, to the measure —Irl used in the derivation of ([T]) and 
retains for general q the various desirable properties of Ikl- Using the more general (H]) leads 
to another nonlinear Schrodinger equation which can be used to calculate leading order energy 
shifts as in Ref.[8j. For unperturbed states with nodes in a symmetric potential the correction 
may be written 
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For integer values of q one may evaluate the function J(?7, q) explicitly but this becomes time 
consuming for large q values. Thus we evaluated it numericaly for a range of real < g < 1. 
The results are as follows. 

For fixed q > 1 and < r/ < 1: J has one local maximum at r]max where it is positive, 
one global minimum at fjmin ^ Vmax where it is negative, and it becomes zero at three points 
= % < ^1 < ^2 < 1- Vmax is between zero and rji while r]m,in is between r/i and 772- J diverges 
to positive infinity as 77 — > 1. The values of 7/i,??2 decrease as q increases. The minimum value 
of J also increases (becomes less negative) as q increases. In the limit g — > 1, 772 — > 1 and thus 
J vanishes at r/ = 1 as we saw in the last section. 

For fixed < g < 1 and < rj < 1: J has one global maximum at rjmax where it is positive, 
one global minimum at Ty^^n ^ Vtriax where it is negative, and it becomes zero at two points 
= % < "^1 < 1- Vmax is between zero and r]i while rjmin is between rji and 1. J is negative at 
r] = 1. The value of r/i decreases as q increases. The minimum value of J decreases (becomes 
more negative) as q increases. 

The results of the last two paragraphs show that J{r], q) reaches its minimum (which is 
negative) for q = 1 and for rj = 0.8; see Fig.(l). Thus using the motivation of the last section, 
minimising the energy of the system, or equivalently minimising our information about the 
system, leads to the q = 1 nonlinear equation ([T]) as the preferred choice, at least when analysing 
static properties of states with definite parity. 

So far we have focussed on states with nodes. For states without nodes the energy shifts are 
0{LP'). For the simple harmonic oscillator we have explicitly 

5E{q) = rf{rf{q + l)(g + 2) - A^{q + 1) + 2)^ + O(L^) . (7) 

This has a global minimum at q ^ 1.1 and 77 ~ 0.76, very close to the values for the states with 
nodes. 
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Figure 1. The magnitude of the minimum value of J{r]-,q) for each q is plotted 



4. Other results and outlook 

Some non-perturbative aspects of the nonlinear Schrodinger equation ([T]) were studied in jT2j . 
In [13] we used Eq.([T]) to study the issue of singularity avoidance in quantum cosmology and 
further work in this direction is in progress. We have also used information theoretic ideas to 
construct integrable nonlinear Schrodinger equations in Ref . [14j . An application of the nonlinear 
Dirac equations of Ref. [5] is in Ref. [E] . 

We welcome correspondence from interested readers. 
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